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We describe the design, performance, sensitivity and results of our recent experiments using the Australia 
Telescope Compact Array (ATCA) for lunar Cherenkov observations with a very wide (600 MHz) bandwidth 
and nanosecond timing, including a limit on an isotropic neutrino flux. We also make a first estimate of the 
effects of small-scale surface roughness on the effective experimental aperture, finding that contrary to expecta- 
tions, such roughness will act to increase the detectability of near-surface events over the neutrino energy-range 
at which our experiment is most sensitive (though distortions to the time-domain pulse profile may make identi- 
fication more difficult). The aim of our "Lunar UHE Neutrino Astrophysics using the Square Kilometer Array" 
(LUNASKA) project is to develop the lunar Cherenkov technique of using terrestrial radio telescope arrays 
for ultra-high energy (UHE) cosmic ray (CR) and neutrino detection, and in particular to prepare for using the 
Square Kilometer Array (SKA) and its path-finders such as the Australian SKA Pathfinder (ASKAP) and the 
Low Frequency Array (LOFAR) for lunar Cherenkov experiments. 
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I. INTRODUCTION 

The origin of the ultra-high energy cosmic rays (UHE 
CR) — protons and atomic nuclei with observed ener- 
gies above 10 18 eV and up to at least 2 x 10 20 eV — 
is obscured due to their deflection and scattering in cos- 
mic magnetic fields. This makes the flux of all but the 
highest energy CR appear almost isotropic with respect 
to the Galaxy regardless of their source, so that measure- 
ments of arrival directions cannot reliably be used for 
source identification. At the highest energies, the de- 
flection is less, and this allows the possibility of 'see- 
ing' nearby UHE CR sources. Arrival directions of 
UHE CR detected by the Pierre Auger experiment above 
5.6 x 10 19 eV are found to be statistically correlated with 
positions of nearby active galactic nuclei (AGN), which 
are in turn representative of the large-scale distribution of 
matter in the local universe [|lj]. However, the flux is ex- 
tremely low, and so the nature of the sources of UHE CR 
within this distribution remains at present unresolved. 

An alternative means of exploring the origin of UHE 
CR is to search for UHE neutrinos. As first noted by 
Greisen [||] and by Zatsepin & Kuzmin [Q], cosmic rays 
of sufficient energy will interact (e.g. via pion photo- 
production) with photons of the 2.725 K cosmic mi- 
crowave background (CMB), with the resulting energy- 
loss producing a cut-off in the spectrum (the 'GZK cut- 
off) at around ~ 10 20 eV from a distant source. These 
same interactions produce neutrinos from the decay of 
unstable secondaries. Several experiments [PJ-EJ have 



reported UHE CR events with energies above 10 20 eV, 
and therefore a flux of these "cosmogenic neutrinos" is 
almost guaranteed. 

Significant information on the CR spectrum at the 
sources is expected to be preserved in the spectrum of as- 
trophysical neutrinos JlC| ] which varies significantly be- 
tween different scenarios of UHE CR production. These 
include acceleration in the giant radio lobes of AGN, the 
decay of super-massive dark matter particles or topologi- 
cal defects, and Z-burst scenarios, the last of which have 
already been ruled out by limits placed on an isotropic 
flux of UHE neutrinos flj"i| |T2[ ], Of course, neutrinos 
are not deflected by magnetic fields, and so should point 
back to where they were produced, with even a single de- 
tection allowing the possibility of identifying the source 
of UHE CR. Here we emphasise that in all models of 
UHE CR origin we expect a flux of UHE neutrinos. See 
refs. [ pj[ |l4| ] for recent reviews of UHE CR production 
scenarios and radio techniques for high-energy cosmic 
ray and neutrino astrophysics. 



A. The lunar Cherenkov technique 

A high-energy particle interacting in a dense medium 
will produce a cascade of secondary particles which de- 
velops an excess negative charge by entrainment of elec- 
trons from the surrounding material and positron anni- 
hilation in flight. The charge excess is roughly propor- 
tional to the number of particles in electromagnetic cas- 
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cades, which in turn is roughly proportional to the en- 
ergy deposited by the cascade. Askaryan [ |l5[ ] first noted 
this effect and predicted coherent Cherenkov emission in 
dense dielectric media at radio and microwave frequen- 
cies where the wavelength is comparable to the dimen- 
sions of the shower. At wavelengths comparable to the 
width of the shower, the coherent emission is in a nar- 
row cone about the Cherenkov angle 8c = cos _1 (l/«) 
(n is the refractive index), while for wavelengths com- 
parable to the shower length, the coherent emission is 
nearly isotropic. The Askaryan effect has now been ex- 
perimentally confirmed in a variety of media |p^-[T^|, 
with measurements of the radiated spectrum agreeing 
with theoretical predictions (e.g. ref. Jl9[). If the inter- 
action medium is transparent to radio waves, the radia- 
tion can readily escape from the medium and be detected 
remotely. Since the power in coherent Cherenkov emis- 
sion is proportional to the square of the charge excess, 
i.e. to the square of the energy deposited, extremely high 
energy showers should be detectable at very large dis- 
tances. 

The lunar Cherenkov technique, first proposed by 
Dagkesamanskii and Zheleznykh [p0||, aims to utilize 
the outer layers of the Moon (nominally the regolith, a 
sandy layer of ejecta covering the Moon to a depth of 
~10 m) as a suitable medium to observe the Askaryan 
effect. Since the regolith is transparent at radio fre- 
quencies, coherent Cherenkov emission from sufficiently 
high-energy particle interactions (specifically, from UHE 
cosmic ray and neutrino interactions) in the regolith 
should be detectable by Earth-based radio-telescopes 
(Askaryan's original idea was to place detectors on the 
lunar surface itself). First attempted by Hankins, Ekers 
& O' Sullivan [^2| ^3p using the Parkes radio telescope, 
the Goldstone Lunar UHE neutrino Experiment (GLUE) 
at the Goldstone Deep Space Communications Complex 
[11], the experiment at Kalyazin [|TJ, and the NuMoon 
[ 24] experiment at the Westerbork Synthesis Radio Tele- 
scope (WSRT) have subsequently placed limits on an 
isotropic flux of UHE neutrinos. Our own project, LU- 
NASKA, aims to develop the lunar Cherenkov technique 
for future use of the SKA. Observations continue at both 
WSRT [ p5| ] and with our own project using the ATCA 
and using the Parkes radio telescope (to be reported else- 
where). The technique has been the subject of several 
theoretical and Monte Carlo studie s Jl^ , p6|~p9[| together 
with our own recent work [23 3(], |31fl. 

Future radio instruments will provide large aperture 
array (AA) tile clusters and arrays of small dishes with 
very broad bandwidths, with both factors allowing very 
strong discrimination against terrestrial radio frequency 
interference (RFI). The culmination of the next gener- 
ation of radio instruments will be the Square Kilome- 
tre Array, to be completed around 2020, with smaller 
path-finders such as ASKAP (Australian SKA Pathfinder 
133]) to be built in the intervening period. In the mean- 



time, we have been performing a series of experiments 
with the Australia Telescope Compact Array [p3|], an ar- 
ray of six 22 m dishes which were undergoing an up- 
grade to an eventual 2 GHz bandwidth at the time of the 
observations described here. Lunar Cherenkov experi- 
ments with these instruments, together with those pro- 
posed for LOFAR p5|], represent the foreseeable future 
of the technique. We emphasize that the lunar Cherenkov 
technique is very different to conventional radio astron- 
omy and requires non-standard hardware and signal pro- 
cessing as it is necessary to detect nanosecond-duration 
lunar Cherenkov radio pulses coming from a region too 
large (the apparent diameter of the Moon is 0.5°) to im- 
age in conventional ways at nanosecond time resolution. 
Such pulses suffer dispersion in the Earth's ionosphere, 
and our experiment is the first to correct for this in real- 
time. 

In the present paper we describe our recent experi- 
ments using the ATCA using the lunar Cherenkov tech- 
nique. We start by giving an overview of the experi- 
ment, the part of the moon targeted and observing times 
which were chosen in order to observe the Galactic Cen- 
ter and Centaurus A (UHE neutrino flux limits for these 
sources are reported in ref. [Q]), the antennas used, spe- 
cialized hardware, triggering and signal processing. We 
then discuss the effects of dead-time, our finite sampling 
rate and our approximate de-dispersion on the detec- 
tion efficiency and effective observation time. Finally, 
we present a new limit to an isotropic UHE neutrino 
flux, and discuss why it is important despite better lim- 
its existing from the Antarctic Impulsive Transient An- 
tenna (ANITA) and the Radio Ice Cherenkov Experi- 
ment (RICE). In appendices, we outline the effects of the 
two major theoretical uncertainties in our aperture and 
limit calculation being the UHE neutrino cross-section, 
and the small-scale lunar surface roughness. In the latter 
case, a new approximate treatment is described, which 
we use as well as the standard calculation methods to 
determine the experimental apertures and limits given in 
the main body of the paper. 



II. DESCRIPTION OF THE EXPERIMENT 

The Australia Telescope Compact Array (ATCA - see 
Fig. |l]) is an aperture synthesis telescope located at lati- 
tude -30° near Narrabri, NSW, Australia. It consists of 
six identical 22 m dishes. One antenna is fixed in po- 
sition, while the other five can be moved along a 3 km 
East- West baseline, and also a ^ 150 m North-South 
baseline. We used three of the moveable dishes on the 
East- West baseline for our observations and used triple 
coincidences with the correct relative timing to identify 
pulses coming from the direction of the Moon 

The ATCA was chosen as an ideal SKA test-bed for 
the lunar Cherenkov technique because: the antennas 



FIG. 1 : A photograph of three of the six ATCA antennas. For the observations described here, the distance between the antennas 
was greater. 



have a size comparable to that expected for the SKA; 
the beam-size matches the lunar disk at 1 GHz; it pro- 
vided us with 600 MHz of bandwidth (to be upgraded to 
2 GHz); because it is an array it provides strong timing 
discrimination against terrestrial RFI; and because it can 
give a large aperture sensitivity while seeing the entire 
moon p5[ ] . Full details of the experiment and observa- 
tions are given by James |36|. 



A. Observations 

The observations described here cover two observing 
periods: February and May 2008, all at the ATCA. Table 
| specifies the array configuration used in each observa- 
tion period. The baselines used were a compromise be- 
tween baselines long enough to resolve correlated ther- 
mal emission from the Moon, and short baselines to re- 
duce the search window in the time domain for pulses 
coming from any location on the Moon. For these ini- 
tial observations we implemented our special hardware 
on only three of the ATCA antennas: CA01, CA03 and 
CA05. 



1. Observation times and antenna pointings 

We had two main considerations in choosing our ob- 
servation times, the first being to confine the observations 
to within the approximate hours of 10 pm to 6 am, in or- 
der to have a stable ionosphere (see Sec. II D 1 ). The re- 



quirement of the Moon being visible meant that this gave 
us a window of perhaps five days in every 29.5-day syn- 
odic lunar month where the Moon would be sufficiently 
visible during this period to warrant observations. The 
second requirement was that the Moon be within ~ 35° 



of particular regions of sky of interest [[31|. This occurs 
once per 27 .3 -day lunar orbit for any given region, so that 
combined, we typically had three good and two marginal 
periods of a few nights each year in which to observe any 
given source. 

The February 2008 run was tailored to 'target' a broad 
(> 20°) region of the sky near the Galactic Center, har- 
bouring the closest super-massive black hole to Earth, 
and a potential accelerator of UHE CR. The Galactic 
Center may also be a source of UHE CR and neutri- 
nos through the decay of massive particles in its dark 
matter halo (see [ |37| ] and references therein). Prelimi- 
nary calculations showed that for beam-sizes similar to 
that of the ATCA, the greatest total effective aperture 
(and hence sensitivity to an isotropic or very broadly- 
distributed flux) is achieved when pointing the antennas 
at the center of the Moon, so all the limb is at approxi- 
mately the half power point of the antenna beam. Since 
any UHE neutrino flux from this region is likely to be 
broadly-distributed, we used this pointing for these runs. 
Our May 2008 observing period targeted Centaurus A 
only, the nearest active galaxy which could potentially 
account for some of the UHE CR events observed by the 
Pierre Auger observatory and to achieve the maxi- 
mum sensitivity to UHE neutrinos from this source we 
pointed towards the portion of the lunar limb closest to 
Centaurus A so it is observed at full sensitivity. 



B. Specialized hardware 

The background signal above which any genuine 
nanosecond-duration lunar Cherenkov pulses had to be 
detected consisted of two components: random noise 
fluctuations, mainly thermal emission from the lunar disk 
and to a lesser extent system noise and Galactic plane 
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TABLE I: Positions of the antennas used during the LUNASKA lunar observations, and the baselines. 



Date 


Configuration 


Antenna stations 


Baselines (m) 






CA01 CA03 CA05 


1-3 1-5 


3-5 


Feb 26-28 2008 


750B 


W98 W113 W148 230 766 


536 


May 18 2008 


750A 


W147 W172 W195 


383 735 


352 


May 19 2008 


EW352 


W102 W109 W125 


107 352 


245 



synchrotron emission; and man-made RFI. Fig. g gives 
a diagram of the hardware and signal path at each an- 
tenna. In order to perform a search for short-duration 
lunar pulses, against a background of thermal noise fluc- 
tuations and RFI, we had to build specialized hardware 
to detect and store candidate events in real time. For 
this we used the digital, field-programmable gate array 
(FPGA) based analog-to-digital converters (ADC) used 
in the Compact Array Broad Band (CABB) upgrade [53], 
each of which could digitize and perform simple logic on 
two data streams at a rate of 2.048 GHz. As well as this, 
we required specialized software running on the control- 
room computers to interface with the CABB analog-to- 
digital converter boards, and a hardware method to cor- 
rect for the dispersive effects of the Earth's ionosphere. 



1. Data channels and signal path 

During each observation period, the received signal - 
split into two orthogonal linear polarizations, which we 
arbitrarily label 'A' and 'B' - was processed by our spe- 
cialized hardware at each of three antennas, providing a 
total of six output data streams. By connecting to a moni- 
tor point at the receivers in each antenna we bypassed the 
normal narrow-band ATCA intermediate frequency (IF) 
system and obtained a nominal bandwidth of 600 MHz 
between 1.2 and 1.8 GHz. For May 2008 a high-pass 
filter was added to remove a strong 0.9 GHz (aliased to 
1 . 1 GHz) narrow -bandwidth RFI from a known transmit- 
ter which was perturbing our detection threshold. 

Each polarized data stream was fed through an ana- 
log de-dispersion filter before being sampled by a 
CABB analog-digital converter board. This operated at 
2.048 Gs/s with 8-bit effective precision. Since both the 
100 Mb/s connection from each antenna to the central 
control room (where all antenna signals could be com- 
bined) and the CABB correlator architecture were inad- 
equate to handle the full raw data rate of 2 x 8 x 2048 
megabit/s per antenna, we had to reduce our data volume 
by triggering independently at each antenna and return- 
ing only blocks of data containing candidate events. The 
signal was copied into both a running buffer and passed 
to a real-time trigger algorithm, and on fulfilling the trig- 
ger conditions a portion of the buffer was returned to the 



control room and recorded. During this recording pro- 
cess, the buffer was unable to respond to further triggers, 
and the experiment was temporarily blind to any events. 
The period of this dead-time depended on the length of 
the buffer to be returned. 



2. Trigger logic and levels 

The trigger algorithm was set up to be a simple thresh- 
old trigger at each antenna - if the square of any single 
sample on either the A or B polarization data streams was 
above a certain value, both polarizations were returned. 
The thresholds were adjusted occasionally to keep the 
trigger rates on each receiver output constant at approx- 
imately 40-50 Hz corresponding to ~ 5a, where we use 
a as short-hand for the rms voltage Vnns in the output 
channel. Even with our 8-bit sampling we were barely 
able to adjust these threshold with sufficient precision 
(< O.lo increments) while maintaining a reasonable dy- 
namic range for any detected event. The gain was ad- 
justed to give an RMS sampler output of approximately 
10 ADC digitization units (a.d.u.) and hence a maximum 
of 12. 8o (128 a.d.u) before saturation for an 8-bit signal. 



C. Dead-time and efficiency 

A certain degree of dead-time loss is suffered for every 
trigger. As all three antennas need to be 'on' to record an 
event, it is important to avoid setting the thresholds too 
low (trigger rates too high) as this can make the effective 
observation time negligible. This dead-time can be eas- 
ily measured by setting the thresholds to zero and record- 
ing the maximum trigger rate for a given buffer length. 
Such measurements were performed at each observation 
period, and the results are recorded in Table We see 
that for a buffer length of 256 samples the dead-time per 
trigger was approximately 1 ms. 

The efficiency of the experiment can be defined as the 
time-fraction when all three antennas are sampling and 
ready to trigger. For a sampling rate r, (Hz) on antenna i, 
maximum rate R{, and purely random trigger events, the 
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FIG. 2: Diagram of the signal path at each antenna. Abbreviations are: 'LNA': low noise amplifier, 'ADC: analogue-digital 
converter. 



efficiency £, is given by: 

5 = n,(i-£) (i) 

where the i multiplies over all three antennas. In Fig. 
H we plot the trigger rates for all antennas for the 18 
May 2008 observations and the efficiency £, calculated 
as in Eq. [IJ assuming a constant Rj = 1040 from Table 
H corresponding to our buffer length of 256 samples. A 
short-duration burst of RFI is evident at UT 12:50 in Fig. 
||, as is a large increase in the background between UT 
15:00 and UT 17:00. During these periods of intense 
RFI the efficiency (upper dashed line) is significantly 
reduced. The effective observation time t e s was deter- 
mined by integrating the efficiency over the observation 
time / bs, and this is given in Table ^together with the 
average efficiency % = ? e ff/? bs- 

D. Ionospheric dispersion 

For our experimental bandwidth of 600 MHz cen- 
tered at 1.5 GHz, the effects of dispersion in the Earth's 
ionosphere are significant. The dispersion is due to 
a frequency-dependent refractive index caused by free 
(ionized) electrons in the ionosphere. Using the stan- 
dard measure for the number of electrons (total electron 



content units, or TECU: 10 16 e /cm 2 ), the time-delay 8f 
relative to a vacuum for a frequency v is given by Eq. (Q): 

8f = 1.34 10~ 7 TECUv~ 2 . (2) 

Of more use is the dispersion At over a bandwidth Av, 
given by 

At = 1.34 10- 7 TECU (v m 2 n -v m 2 x ) (3) 
At « 2.68 10~ 7 TECUAvv~ 3 (4) 

if Av <C V, where v is the mean frequency (v m ; n + 
V m ax)/2. Note that in terms of phase delay, the correction 
goes as v~ 2 . 

1. De-dispersion filters 

Implementing a digital de-dispersion filter running at 
this speed was too difficult at the time, so we used ana- 
log de-dispersion filters designed by Roberts j4C|]. Each 
filter was a variable-width waveguide of approximately 
one meter in length constructed as a spiral for compact- 
ness (Fig. ^), with the output being the continuous sum 
of reflections along the length. Upon reflection, high fre- 
quencies experienced a greater delay than low frequen- 
cies, with the design such that this cancelled out the 
delay due to ionospheric dispersion at low frequencies. 
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TABLE II: Maximum trigger rates (Hz) as a function of buffer length for 18 May 2008. 



Buffer length 16256 


8192 4096 2048 


1028 


512 


256 


128 64 


Trigger rate 22 


42 83 163.5 


317.5 


581 


1040 


1690 2450 



500 
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100 
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80 



60 



40 



20 







12 13 14 15 16 17 
Hours UT (May 18th) 



18 



FIG. 3: (color online) Trigger rates for the three antennas (three lower curves) and efficiency (upper curve) from Eq. (|lj) for the 
night of May 18*^2008. The increased trigger rate at UT 15:00-17:00 results from a 'Type 1' noise feature caused by an unknown 
source - see Sec. [v| and the top-centre of Fig. |Tc| . 



TABLE III: Raw observation time f (, s (minutes), mean effi- 
ciency \ (%), and effective observation time f e g (minutes), for 
all observation periods. 



Date 




\ 


r eff 


26 Feb. 2008 


239 


86 


204 


27 Feb. 2008 


319 


87 


277 


28 Feb. 2008 


314 


87 


274 


17 May 2008 


324 


69 


224 


18 May 2008 


376 


73 


274 


19 May 2008 


440 


72 


316 



Thus an in-phase signal (e.g. coherent Cherenkov radi- 
ation from a UHE particle interaction in the Moon) en- 
tering the top of the ionosphere should appear in-phase 
after de-dispersion, provided the correct dispersion mea- 
sure was used. 

We used NASA data [ f4l| ] to predict typical values 
of the ionospheric dispersion at Narrabri, with results 
reported in Ref. [p9|]. Since we were still near solar 
minimum these results showed that the ionosphere over 
Narrabri was comparatively stable between the hours 
of 10 pm and 6 am, with low and predictable verti- 
cal total electron content (VTEC) measures around 7 ± 
1.3 TECU. This corresponds to a differential vertical de- 
lay of 3.6 ns over a 1.2-1.8 GHz bandwidth. Since the 
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(a) 




FIG. 4: Printed circuit board layout of analog de-dispersion 
filter. Note the variation in the width of the waveguide along its 
length. The physical size is approximately 25x25 cm. 



actual delay will depend on the slant angle, we chose 
to build the filters assuming a 5 ns delay over the band, 
i.e. a TEC along the line of sight (slant TEC, or STEC) 
of 10 TECU; this is also equivalent to the mean VTEC 
of 7 TECU with a lunar elevation of 47°. Therefore we 
expected to lose some sensitivity when the Moon was di- 
rectly overhead, and also very near the horizon. The sen- 
sitivity lost due to deviations of the actual VTEC from 
the mean and variations in lunar elevation is discussed 
in Sec. |ni| Fig. ||(a) shows the expected effect of de- 
dispersion on pulses of different origin: satellite bounce 
(of RFI), solar system (lunar Cherenkov) and terrestrial 
(RFI). These have been modeled respectively as a flat fre- 
quency spectrum dispersed once (dispersed twice then 
de-dispersed once), a v 2 spectrum expected for coher- 
ent Cherenkov emission unchanged (dispersed once then 
de-dispersed once), and a flat frequency spectrum de- 
dispersed once. In all cases the Fourier inverse has been 
taken over 1.2-1.8 GHz, i.e. the pulse has been band- 
limited to 1.2-1.8 GHz. 
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FIG. 5: (a) Predictions of an observed impulse of different ori- 
gins, (b) A typical noise-calibration pulse - in this case, from 
CA01 A on Feb 26 th . The small pre-pulse at —35 ns is caused 
by a spurious reflection from the filter connection point. 



switching of the noise-calibration diode. Note that, as 
expected, the pulse in Fig. ||(b) closely resembles that 
expected for a terrestrial source. 

Since the noise diode could be turned on and off from 
the control room, it could be used to generate approxi- 
mately coincident false triggers between the three anten- 
nas. The time difference in the noise diode switching be- 
tween antennas was small (< 1 /js) and the scatter about 
this offset was ~ 200 ns. Using this simple procedure to 
generate approximately coincident false triggers became 
a useful part of our experimental procedure as discussed 
in Section IV-A. 



E. The noise diode 

A noise-calibration diode is located in the receiver 
and used to calibrate the power through measuring the 
system temperature, r sys , during normal ATCA observa- 
tions. During testing, we observed unexpected strong 
pulses at a rate of approximately 8 Hz (Fig. ||(b)). It 
was discovered that these pulses were generated by the 



III. SAMPLING AND DISPERSIVE EFFECTS ON 
SENSITIVITY 

The effects of loss of sensitivity due to a finite sam- 
pling rate (compared to an infinite sampling rate) and de- 
dispersion depend on the frequency spectrum of the radio 
pulse and the bandwidth of the detector. The expected lu- 
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nar Cherenkov pulse spectrum depends on many factors, 
including lunar surface roughness and orientation, the di- 
mensions and direction of the electromagnetic cascade in 
the lunar regolith, and on the neutrino energy - these are 
discussed in detail elsewhere [ 30 1 . The range of possible 



spectra is very broad. However, for our purposes we con- 
sider two extreme cases. We are concerned with the rel- 
ative strengths of the high-frequency and low-frequency 
components. Near the minimum detectable cascade en- 
ergy for our experiment (10 20 eV - see Sec. |l5]), only 
fully coherent emission will be detectable. The electric 
field spectrum (V/m/MHz) will therefore have the form 
E(v) — Ay, which gives the greatest possible weight to 
the high-frequency component. 

The other extreme is given by a high-energy shower, 
at shallow depth and viewed at a large angle away from 
the Cherenkov angle where the emission from the cas- 
cade is becoming incoherent. In this case we expect 
^(v) = Bvexp(— Cv 2 ) based on approximate fits to the 
spectrum far from the Cherenkov angle for hadronic 
showers For any observed pulse of given total 

power we do not know the shape of the spectrum to ex- 
pect if the pulse is indeed of lunar origin, and so we 
consider the two extreme possibilities. For our 'high 
energy' we take 10 23 eV, since above this range strong 
limits on a neutrino flux from the NuMoon [ ^4| and 
FORTE p0| ] experiments made a detection extremely un- 
likely. For the same reasons, 10 23 eV is also the most 
energetic neutrinos to which we simulate our effective 
apertures and limit, which can at most produce cascades 
of energy 10 23 eV. Setting the peak power in the band- 
width to be the same for the two extreme cases, corre- 
sponding to 10 20 eV and 10 23 eV neutrinos, and taking 
B/A = (10 23 /10 20 ) allows the constant C to be found as 
follows. Since the electric field of a pulse (wave packet) 
at the antenna may be written as 



E(t) 



£(v)e 27C ' v( '- ; °W, 



(5) 



and since the power is proportional to ^(f)! 2 , the peak 
power occurs at t — fo and so is proportional to 



|£('o)l 2 



E(v)dv 



(6) 



Hence, for equal peak powers for the two extreme forms 
of possible spectra we may solve 



„V 2 rV 2 

A vdv = B vexp(-Cv 2 )dv 



(7) 



: 1.2 GHz, and v 2 = 1.8 GHz we 
2 . We shall use these results when 



For B/A = 10 3 , Vi 
find C = 3.515 GHz 
finding the uncertainty in sensitivity due to the unknown 
spectral shape. 
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FIG. 6: (color online) Maximum sampled value as a fraction 
of the true pulse height as a function of the arbitrary phase 
offset for sampling rate of 2.048 Gs/s and two extreme lunar 
Cherenkov pulse types. 



A. The effects of a finite sampling rate only 



Our sampling rate of 2.048 Gs/s was greater than the 
Nyquist rate of 1 .2 Gs/s for our nominal 600 MHz band- 
width, and allowed for perfect reconstruction of the sig- 
nal in the frequency range 1.024 to 2.048 GHz at arbi- 
trary time resolution (assuming no signals outside this 
range). However, experiments (such as this one at the 
ATCA) which are pushing the current FPGA limits for 
high-speed signal processing can only use simple algo- 
rithms based on the pulse height of a single sample as 
captured by the sampling threshold. For any finite sam- 
pling rate (including the Nyquist rate), there will be a 
random offset in the phase of the ADC digitization times 
between the actual peak of the pulse and the sampling 
times. Fig. || plots the peak pulse height as a func- 
tion of the arbitrary phase offset for the sampling rate 
of 2.048 Gs/s as used in our experiment and the two ex- 
treme spectra of lunar Cherenkov emission. The peak 
sampled value is seen to vary by ^30% in the case where 
the Cherenkov emission is fully coherent, and by ~15% 
when the emission is becoming incoherent. Hence, with 
our simple trigger logic, there will be a triggering inef- 
ficiency due to our finite sampling rate causing the peak 
sampled voltage occasionally to be less than the trigger 
level, even if the actual peak voltage in a pulse was above 
it. In future experiments, we will use a more complex 
trigger algorithm which works off multiple sampled val- 
ues in order to reduce this loss, which is a more efficient 
remedy than increasing the sampling rate. 
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B. Effects due to de-dispersion 

In order to quantify the effect of using the constant 
TEC value built into our de-dispersion filters, we take 
pulses of the two extreme types discussed above, dis- 
perse them for different line-of-sight TEC values, de- 
disperse them using the constant TEC value built into 
our de-dispersion filters, and finally simulate sampling 
by the ADC. All possible offsets of the pulse-peak ar- 
rival time with respect to the sampling times, or "base 
phase offsets", were modeled in this process. For each 
combination of intrinsic spectrum, base phase offset, and 
dispersion measure, we calculate the peak pulse strength 
in the time domain. Averaging this over all base phase 
offsets (which will be random) and dividing by the mag- 
nitude of the peak undispersed pulse at zero phase offset 
we obtain the peak signal strength as a function of line- 
of-sight TEC shown in Fig. ^. As well as the sampling 
rate of 2048 GHz used, we also show results for sam- 
pling rates a factor of 2 higher and lower. The upper 
and lower sets of lines are for pulses due to incoherent 
and coherent Cherenkov emission by the lunar cascade 
as a whole, respectively. The pulse for fully-coherent 
Cherenkov radiation is most adversely affected by iono- 
spheric dispersion because the signal is spread over the 
largest frequency range. The rapid oscillations in average 
peak amplitude with changing line-of-sight TEC for the 
coherent Cherenkov pulse is due to the combined effect 
of the sharp band edges and the de-dispersion function. 
The mean values from Fig. || correspond to zero line- 
of-sight TEC for the case of a 2048 GHz sampling rate, 
allowing a comparison between the effects of our finite 
sampling rate and dispersion. We note that the loss in 
sensitivity due to errors caused by using a constant TEC 
value for de-dispersion (typically differing from the true 
value by less than 4 TECU) are less than the finite sam- 
pling rate errors for our observations using a simple trig- 
ger algorithm. 



C. Loss of experimental sensitivity 

The loss of sensitivity due to sampling and disper- 
sive effects can be calculated using the measured values 
of total electron content (TEC) as a function of eleva- 
tion, which were determined after the observations J4l|]. 
The results are given in Fig. |[ Using a linear interpola- 
tion between these points and the known lunar elevation 
gives the line-of-sight TEC (slant TEC, or STEC) mea- 
sure (dotted lines). At low elevations, the line of sight 
will probe a large horizontal distance, so using a con- 
stant VTEC measure may not be appropriate. However, 
since the TEC goes as 1 / sin (elevation) and consequently 
blows up at low elevations, the sensitivity in this regime 
will be low in any case. Combined with the mean losses 
for the two spectra in Fig. ^ the range of losses for the 




20 30 
STEC (TECU) 

FIG. 7: Greatest detected amplitude of the time-domain sig- 
nal relative to the maximum undispersed amplitude, for two 
extreme lunar Cherenkov pulse types and four sampling rates, 
as a function of the difference between the line-of-sight total 
electron content, in TECU, and the constant value used in the 
de-dispersion. 



TABLE IV: Estimated nightly average fractions of the peak sig- 
nal detected (%) for the observations periods indicated. The 
best case corresponds to incoherent Cherenkov emission from 
lunar cascades, the worst case to completely coherent lunar 
Cherenkov emission signals, as discussed in text. 



Period 


February 2008 




May 2008 


Date 


26 th 27 th 28 th Mean 


17 th 


18 th 19 th Mean 


Best (%) 


92 91 91 91 


91 


91 91 91 


Worst (%) 


84 80 81 82 


82 


80 81 81 


Mean (%) 


88 86 86 87 


87 


86 86 86 



experimental periods is calculated as the shaded regions 
in Fig. |8] which are bounded by the extreme spectra of 
coherent/incoherent lunar Cherenkov emission by lunar 
cascades. 

The effect of using a fixed STEC value for the de- 
dispersion is estimated by taking the mean loss over both 
observation time and spectral type, giving equal weight- 
ing to both the coherent and incoherent spectra. The re- 
sulting mean detected signal fractions are given in Table 
IV 



IV. RELATIVE TIMING CALIBRATION WITH 
ASTRONOMICAL POINT SOURCES AND RFI 
SOURCES 

We had counters recording the number of samples at 
each antenna, the value of which was returned with each 
triggered event. These could be converted to clocks accu- 
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FIG. 8: Ionospheric influence over our May 2008 observations. Dashed curves: measured VTEC (TECU) from Ref. Dotted 
curves: STEC (slant TEC along the line of sight). Solid curves: mean fractional loss for coherent (lower curve) and incoherent 
(upper curve) lunar Cherenkov pulses due to ionospheric effects and triggering inefficiency due to our non-infinite sampling rate 
(see Fig. ^): the shading gives the range. Vertical lines indicate the times of the observations, horizontal lines give the fractional 
loss averaged over each observation and over the two extreme cases for the expected spectrum. 



rate to ~ 0.5 ns. However, at the time of our experiment 
these clocks had unknown timing offsets between them, 
which had to be determined to allow a sufficiently rigor- 
ous pulse search. To calibrate the times, we required both 
a common signal in each antenna with some known time- 
delay, and a method to trigger the buffers with sufficient 
simultaneity that enough of the common signal seen by 
all three antennas would be captured to produce a signif- 
icant correlation. We used the noise-calibration pulses 



(see Sec. |IIE] ) as our trigger, and the very bright discrete 
radio sources 3C273 and 3C274 (M87) as our correlated 
signal. 

The QSO 3C273 was chosen since it was the bright- 
est (47 Jy at 1 .4 GHz) point-like source near the Moon 
at the time of the observations, and would thus give a 
strong correlation over all baselines. The radio galaxy 



M87 is brighter (215 Jy at 1.4 GHz) [|33| but resolved on 
our long baselines, and was chosen to maximise the cor- 
related signal over our shortest (CA01-CA03) baseline 
only. 



A. Timing observations 

To calibrate the timing, we pointed the antennas at 
either 3C273 or M87, set the buffer length to maxi- 
mum, switched the noise-calibration on, and set the trig- 
ger thresholds such that we were triggering only off the 



noise-calibration pulses, at roughly 8 Hz. We typically 
observed in this calibration mode for a few minutes at a 
time and thus took of order 2000 pulses, repeating this 
procedure a few times each night. 



The timing offsets, At, measured in samples, are given 
in Fig. H for the February observation periods. The ver- 
tical axis shows the absolute time offsets between anten- 
nas /' and i, A/P., relative to the first calibration after each 
clock reset - i.e. the first data points have been adjusted 
to 0. This adjustment is the timing calibration offset. 
The expectation was that all data would therefore have 
y-values near 0. Obviously this is not the case, and we 
see the time offsets Atfi jump around in multiples of 192 
samples (93.75 ns). 



We eventually discovered that these 192-sample off- 
sets were a hardware fault triggered whenever we 
changed the buffer lengths. In the majority of cases we 
were able to use our log-book and RFI in the data to suc- 
cessfully identify when these offsets occurred and make 
the appropriate correction. In a small number of time 
blocks we still had ambiguity so we searched for event 
coincidences using all possible 192-sample offsets. Be- 
cause the triple coincidence requirement is so strong this 
did not have an impact on the final sensitivity. 
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FIG. 9: Consistency check of correlation times in the February 
data: alignment relative to the first calibration of the February 
period, in (x-axis) chronological order. 



TABLE V: Number of two-fold and three-fold coincidences, 
within an 8000 sample (~ 4 /js) window for each night (the 
count of two-fold events also includes the three-fold events). 



Date 


CA01/03 


C AO 1/05 


CA03/05 


CA01/03/05 


February 26 


6445 


1286 


1533 


449 


February 27 


68894 


39898 


43224 


30051 


February 28 


23296 


8590 


11072 


6781 


May 17 


2344 


1925 


1994 


96 


May 18 


21774 


19445 


20635 


3437 


May 19 


114383 


74311 


71313 


57493 



V. RFI 

After correcting the timing for the 192-sample jumps, 
we looked at coincident triggers within the physically 
possible time range for signals coming from outside the 
array. The number of two-fold and three-fold coinci- 
dences for each observation night are given in Table |V[ 
The most obvious result is the extremely large number 
of two-fold coincidences, and the large number of three- 
fold coincidences compared to two-fold coincidences. 

The expected rate (Hz) of two-fold coincidences, Ry, 
and three-fold coincidences, R135, from purely random 
arrival times is given by 



Rij = RiRjWt 
R l35 = R,RjR k W, 2 
R l35 /Rij = R k W t . 



(8) 
(9) 
(10) 



where Ri is the rate (Hz) of single triggers in antenna i 
and Wt is the time window (seconds) required for a coin- 
cidence. Hence, the ratio between three-fold and two- 
fold coincidences increases with the trigger rate. For 



a maximum trigger rate of 3 kHz and time window 
W t = ±3.906 1(T 6 s (8000 samples), the two-fold rate is 
70 Hz. However, the three-fold rate is only 1.6 Hz, i.e. 
only ~ 2% of the two-fold rate. While at times the ratio 
of two-fold to three-fold coincidences matched this ex- 
pectation exactly, there were also some periods within 
each night where up to 90% of coincidences between 
CA01 and CA05 were also coincident with a CA03 
event, indicating an RFI source. The obvious conclusion 
therefore is that the vast majority of observed three-fold 
coincidences do not occur purely randomly, but rather 
are triggered from a common event with significant time- 
structure. By extension, there will be many such events 
seen only in two antennas, and the same must therefore 
apply to the two-fold coincidences, of which there are 
(generally) many more. 

In Fig. [l(], we have plotted all two-fold coincidences 
for each of the three pairs of antennas over the entire ob- 
serving period for both runs - over 5 x 10 5 data points. 
For each two-fold coincident event the angle, 9, between 
the East direction and the direction of propagation of a 
plane wave fitting the arrival times at the two relevant an- 
tennas is calculated, and cos 8 is plotted (dots) against the 
time of the occurrence of the event (top) February, and 
(bottom) May 2008. Time increases continuously except 
in breaks between days as indicated by the thick dashed 
vertical lines. Note that regions of high 'dot-density' ap- 
pear blue only because the blue dots were plotted last. 
The lunar direction is indicated by the grey dotted line, 
and positions of candidate events (crosses) are marked. 
For May 17th, a network error necessitated many adjust- 
ments to the buffer size, resulting in multiple unknown 
timing offsets to occur during that night. Therefore the 
timing criteria for candidate lunar events have been re- 
laxed to ±5 sample offsets (i.e. ±0.47 /js) on all antennas 
for May 17th, explaining why the four candidate events 
for that night do not lie exactly on the Moon's trajec- 
tory. For May 18 the trial offsets were ±1 sample offsets 
(i.e. ±94 nanoseconds) on all antennas, the system was 
more stable and had two calibrations in agreement, and 
on May 19 (and all of February) the times were com- 
pletely aligned. 

A stationary source of pulse-like RFI producing trig- 
gers over a long period of time will show up as a horizon- 
tal line, while a very brief period of strong, narrow-band 
RFI will be observed as a very large number of coinci- 
dences over a small time range but large vertical extent 
(since the times will be random), i.e. a vertical line. 

Fig. [T(] provides an amazing amount of information. 
The features can be approximately classed as below: 

1 . Short time periods exhibiting a high rate of coin- 
cident triggers for all cos 9, mostly in May (e.g. 
15:35 UT May 17 th , 12:50 UT May 18 th , 10:00- 
17:20 UT May 19 th ). These appear as vertical fea- 
tures which extend uniformly over the full range 
of offsets. 
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FIG. 10: (color online) For each two-fold coincident event the angle, 6, between the East direction and the direction of propagation 
of a plane wave fitting the arrival times at the two relevant antennas is calculated, and cos 9 is plotted (dots) against the time of 
the occurrence of each event for February (top), and May (bottom). Red dots represent coincidences between antennas CA01 and 
CA03, green dots for CA01 and CA05, and blue dots for CA03 and CA05. In dot-dense regions, only blue points show, since they 
are plotted last. The lunar direction (grey dotted line) is also plotted, and positions of candidate events (crosses) are marked. 



13 



2. Short time periods exhibiting a high rate of coin- 
cident triggers for a broad range of cos 8, mostly 
in February (e.g. 14:25 UT February 26 th , 14:40 
and 14:50 UT February 28 th , and 19:40 UT May 
19 th ). These appear as vertical features which do 
not extend over the full range of offsets. 

3. Purely horizontal, typically thin features occurring 
at a characteristic offset, sometimes over many 
days (e.g. the line at cos 9 = -0.88 for 17:00 UT 
on February 26th and cos 9 = 0.92 for 14:30 UT 
on May 19th. 

4. Sloped, typically broad features occurring most 
obviously around 17:00-19:00 UT on February 
27 th , but also on February 26 th . 

Type 1 features are exactly what would be expected from 
a high random trigger rate, with triggers evenly spread 
in time-offsets. The cause of the increased trigger rate 
must be a lowering of the effective threshold by an in- 
crease in the background containing no timing informa- 
tion, probably from narrow-band RFI - an increase due 
to ground temperature or the galactic background would 
be unlikely to produce such short-duration bursts. As 
would be expected, the triples rate is much lower than 
the doubles rate during these times, due to the random 
nature of the trigger times. 

Type 2 features show both a small increase in triggers 
at all offsets corresponding to a random component from 
an increased background, and a large increase in triggers 
spread about a specific offset. Some bands likely corre- 
spond to the nearby towns of Narrabri (East) and Wee 
Waa (West-Northwest) which present potentially large, 
extended sources of RFI. While these may be too weak to 
be detected under normal conditions, a decrease in the ef- 
fective threshold due to the presence of narrow -band RFI 
could make such sources detectable. Another explana- 
tion is a single source of RFI with a broad time-structure, 
with a high likelihood of different antennas triggering 
off different parts of the signal, thereby adding a random 
component with a preferred direction. These features ap- 
pear strongly for both doubles and triples, indicating the 
non-randomness of the timing. 

Type 3 features act exactly as fixed sources of short- 
duration RFI. The fixed geometrical delay results in a 
common time-offset regardless of antenna pointing posi- 
tion (and hence time). 

Type 4 features remain largely unexplained. The locus 
of the coincident triggers (lines with constant slope) are 
consistent with some RFI source moving with constant 
speed in cos 9 (9 the angle w.r.t. the baseline), which is 
strong enough to give a high rate of triple-coincidences, 
and in the far-field, since the delays per unit baseline 
match. Since multiple features are seen at once, there 
must be either many such sources all moving in uni- 
son, or many multiple reflections keeping the same (and 



extremely large) angular offsets over a broad period of 
time. Also, the apparent motion is at the sidereal rate, 
but in the opposite direction. One suggestion is that a 
far-field RFI source is being observed over multiple sig- 
nal paths due to tropospheric ducting which is associated 
with inversion layers. Reflections off the antennas them- 
selves cannot explain the rate of change of delay being 
proportional to the baseline length, nor is the antenna size 
of 22 m sufficient to produce more than a ~ 70 ns change 
in delays. We can be certain however it is not an equip- 
ment fault due to the presence of the aforementioned type 
3 feature during this time period. 

The apparent anti-sidereal motion of the features may 
result from either the real motion of far-field sources of 
RFI, or a series of reflections off an extended object, al- 
lowing each reflection point to move smoothly with time. 
In the former case, a possible candidate is a set of satel- 
lites in a medium Earth orbit (altitude ~ 20,000 km), 
which should move west-to-east across the sky at ap- 
proximately the correct rate. This orbital altitude is occu- 
pied primarily by navigational satellites, such as those of 
the Global Positioning System (GPS), with 'LI' and 'L2' 
carrier frequencies of 1575.42 and 1227.6 MHz respec- 
tively. The positions of the GPS satellites over the period 
of the experiment were checked from public ephemeris 
data [43 1, and found to exhibit the expected anti-sidereal 
motion, but they did not match the positions of the fea- 
tures. A more extensive search of all satellite positions 
might uncover a suitable candidate however. 



A. Source identification with three-fold triggers 

If any on-site RFI sources are found, they could be 
deactivated and/or shielded in time for follow-up exper- 
iments, since in many cases these events dominate our 
trigger rates and limit sensitivity. Given three antenna 
positions on an East-West baseline, we can solve for 
the source position to within a North-South ambiguity, 
since an event some distance North of the baseline would 
produce exactly the same time structure if its location 
were directly South of the baseline by the same distance. 
We break the three-fold coincidences into two types of 
events: near-field and far-field. 

Using the timing offsets, a search for both far-field 
and near-field events was performed for each block of 
data in each of the February and May observation pe- 
riods. The majority of near-field solutions occur in the 
very near-field (within 1 km of the antennas), and the 
rates of both near-field and far-field events are highly 
variable. In many cases point-like sources of RFI are 
seen, both in the near-field and far-field, and it makes an 
interesting game trying to align the positions of possi- 
ble RFI sources with those detected. Probable sources of 
RFI that we identified in this way include the residence, 
the control-building/lab, the solar observatory, the lodge, 
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and either or both of the Ionospheric Prediction Service 
center or visitors center, all of which are on the ATNF 
site at Narrabri. 



VI. SENSITIVITY CALIBRATION 

To simulate the sensitivity of this experiment and place 
limits on a flux of UHE neutrinos, trigger levels at each 
antenna in terms of real quantities must be known. The 
usual specification for an experiment such as this is the 
detection threshold in V/m/MHz (a V/m threshold 
divided by the bandwidth) in a given polarization just 
prior to being received by the antenna. For instance, the 
GLUE threshold was a maximum field strength per band- 
width of £thr = 1-23 1(T 8 V/m/MHz, calculated by tak- 
ing the threshold of 6.46 10~ 9 V/m/MHz in each circu- 
larly polarized data channel, and accounting for vacuum- 
receiver transmission and the splitting of power between 
polarizations [p4|], To calculate our V/m/MHz thresh- 
old, we first had to perform our own calibration of the 
antenna gain as a function of frequency (bandpass cal- 
ibration), since the automated ATCA measurement of 
r sys using the injected noise source only applies to the 
standard signal path and over a small frequency range. 
Hence, another method had to be used, as described in 
the section below. Also discussed below are the effects of 
ionospheric dispersion, which while approximately cor- 
rected for, still reduced our sensitivity to some degree. 



A. The calibration function k(v) 

For a wide-band experiment such as this, the signal 
is expected to change significantly in strength over the 
band. Note that unless otherwise stated, Av refers to 
the total bandwidth from 1.024-2.048 GHz recorded, al- 
though the sensitivity outside the range 1 .2-1 .8 GHz will 
be minimal. Our sensitivity will change as a function 
of frequency due to the antenna, the receiver amplifier 
bandpass, and the de-dispersion filter. Given that our 
threshold is set in terms of our 8-bit sampling, we need 
to be able to convert from a signal at the antenna re- 
sulting from coherent Cherenkov emission from a lunar 
cascade, E(y) (V/m/MHz), to the value of the received 
buffer b(t) at the peak of the pulse (t = f pea k)- From 
here on, we simply call the values of b "ADC digitiza- 
tion units", or "a.d.u.", and all frequencies are in units 
of MHz. The relationship between E(y) and £>(?p ea k) 
involves an unknown (but determinable) function k(v), 
which is required to calculate the peak signal height (ar- 
bitrarily, occurring at time t = 0), defined as below: 



Mfpeak) (a.d.U.) 



k(v)E(v)dv. (11) 



This gives the conversion between real field strength at 
the antennas and the measured units in the CABB ADC 
boards. For coherent pulses away from the peak, and for 
incoherent signals at all times, the integral on the RHS 
of Eq. |ll] should include a phase factor e 2n,vl , i.e. it is a 
Fourier transform. Therefore k(v) can be more simply 
defined with respect to the Fourier transform b(v) of bit ) 
as per Eq. O: 



b(v) = Jt(v)£(v). 



(12) 



For simplicity, a more useful measure is k, being the 
mean over the bandwidth Av between v m j n and V max : 



1 /'Vmax 

k = — k(v)dv. 
Avi Vmin 



(13) 



The sensitivity of the experiment E^, defined in terms 
of a threshold electric field strength per unit bandwidth 
(for the simplest case of a flat-spectrum pulse), can then 
be calculated by knowing the trigger threshold b^ using 
Eq.pl: 



£thr(v) (V/m/MHz) 



?thr 



(0 



kAv 



(14) 



for Av in MHz. Therefore, in this section we calculate 
separately k(v), k, and hence for each data channel 
(antenna and polarization) over the entire observation pe- 
riod. 

In order to calculate k(v), a measurement of a known 
flux F(v) (W/m 2 /Hz) is required. For an incoherent sig- 
nal (random phases), the relationship between F(v) and 
the electric field over a given bandwidth is given by: 



F(v)dv = E^/Zo, 



(15) 



where Zq = fj^c is the impedance of free space, and the 
RMS electric field is measured by the antenna system as: 



1 

At 



E{t) 2 dt 



(16) 



Parseval's theorem for the E(t) ■<-» E(v) transform tells 
us that: 



E{tfdt 



E(yfdv. 



(17) 



Substituting Eq. dl7| ) into (|T6|), and thence into Eq. (Jl5|), 
we arrive at Eq. (tt8fi: 



•'v„ 



F(v)dv 



Z Q At 



E(v) 2 dv. (18) 



Since this relationship holds for an arbitrary bandwidth, 
the integration can be eliminated (e.g. let v max — > V m ; n ), 
giving: 



F(V) 



1 

ZoAt 



m 2 



(19) 
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Using Eqs. (11) and (p^), the flux F(v) as seen by the 
relevant data channel can then be related to the required 
calibration constant k(y): 



F(y) 



1 



b(y) 



Z At \k(v) 
which can be rearranged to give k(v): 



*(v) 



b(v) 



v /F(v)Z Af 



(20) 



(21) 



Eq. ( pi] ) states that we can obtain k(v) from a known flux 
F(v) and the Fourier transform b(v) of the correspond- 
ing sampled output. Note that the output is actually sam- 
pled discretely, so b(v) is obtained indirectly, through a 
discrete Fourier transform. Note also that for incoher- 
ent thermal emission, we do not need to track the phase 
change in each antenna over the bandwidth, i.e. we are 
interested only in the magnitude of k(v). 



B. The Moon as a flux calibrator 

For our calibrator, we chose the Moon. The lunar tem- 
perature Tm is stable to within a few degrees over the lu- 
nar cycle at approximately 225 K in the 1-2 GHz range 
(see Ref. |p5]]). There are small errors in this assumption 
due to the variation (1-2%) with lunar cycle, comparable 
variation across the band, and variations across the disk 
of the moon, and small polarization effects. Combining 
these errors, this method should be accurate to within 5% 
or better, which is acceptable. 

Under these approximations, the lunar flux Fm(v) (Jy) 
captured by the beam (it will be half this in any given 
polarization channel) is given by: 



F M (y) 



2k h T M v 2 



<B[Q(Cl,Cl p ),v]dQ. (22) 



where kt, IS the Boltzmann constant, and $(8,v) is the 
beam-power pattern of ATCA telescope dishes [ p^| at 
frequency v at angle 8 to the telescope pointing direc- 
tion fo p , (i.e. 8 is the angle between directions fo and 
fop), and Q.m is the solid angle subtended by the Moon. 
The beam-power pattern of ATCA deviates only slightly 
from the Airy pattern of a 22 m diameter aperture. In the 
present observations, fo p was either the direction towards 
the center of the moon or towards the lunar limb. 



C. Measurements 

To obtain b(v) by observing lunar thermal emission 
we took an unbiased sample of data pointing both on 
and off the Moon by setting the trigger level to zero, i.e. 



maximally triggering. The received flux Fm{v) from the 
Moon can be detected by subtracting the measured band- 
pass b g(v) when pointing away from the Moon from 
the bandpass b on (v) when pointing at the Moon's cen- 
ter. The pointing-position for the off-Moon data was a 
position at similar galactic latitude far from any strong 
sources in the ATCA catalog. We set the buffer lengths 
to maximum for this procedure, since then the product 
of trigger rate and buffer length is largest, and also we 
obtain the best spectral resolution. This was done once 
every time the configuration was changed. 

Each of the Nt recorded buffers was discrete-Fourier- 
transformed to produce b(y) (typically Nt ~ 5000). The 
resulting spectra are squared and then averaged over all 
the buffers recorded for each calibration period/target 
taken. Each averaged spectrum is then cleaned with a 
very simple cleaning algorithm to remove the worst of 
the RFI, which simply sets the power of all RFI spikes 
above a running threshold to zero, and the subsequent 
analysis ignores them. An example of the raw and 
cleaned spectra is given in Fig. |ll} 

The squaring, summing, and cleaning process was 
repeated for both the off-Moon and on-Moon (center 
and limb) spectra for each antenna/polarization, and the 
off-Moon power-spectrum is subtracted from the cor- 
responding on-moon spectra. Taking the square root 
gave the required \b(v) \ corresponding to the lunar con- 
tribution as required for Eq. (|2l|). Note that we make 
the approximation \b(v)\ ~ b(v) assuming that the de- 
dispersion filter keeps the phase constant across the band. 
This was then divided by the product ^/ZoFm{v)At to 
give |*(v)|. 



D. Results of the calibration 

1. Fitting for k(v ) 

In order to characterize k(v) in a meaningful way, a 
piece-wise linear approximation to k(y) was performed. 
Fig. |l| shows the fits for CA01 A in February 2008 - 
four fits have been used, with different frequency ranges 
for each of the February and May periods. Also shown 
is the mean k(v), which has been fitted to the bandwidth 
1.1-1.8 GHz in the case of February. For simulation pur- 
poses, the piecewise-linear fit to the bandwidth was used. 
Note that we have smoothed over the oscillations in k(v), 
which are caused by interference between the dispersion- 
corrected signal from the filters, and a small reflection 
from the filter connection point - this reflection is also 
the cause of the pre-pulse observed in Fig. 0. 

While some sensitivity is not included by limiting the 
range of the fitted bandwidth, including this range in 
the fit would artificially reduce k(v) at lower frequencies 
where a signal is more likely to be observed. Conversely, 
taking the fit below the low-frequency cut-off would have 
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FIG. 1 1 : (color online) Raw (top) and cleaned (bottom) spectra for the off-Moon pointing of May 18 on CA05. 



led to an overestimate of the sensitivity at low frequen- 
cies where the signal is stronger. Unsurprisingly, the 
ranges which gave a good trade-off between these effects 
and artificially reducing the effective bandwidth Av were 
close to the nominal bandwidth of 1 .2-1 .8 GHz. 

Only one calibration was performed in February, and 
no cross-checking was possible. In May, it was found 
that while fits from both May 19 th measures and the limb- 
off-Moon fit of May 18 th were always in good agree- 
ment (±3%), the lunar center/off-Moon fits of May 18 th 
were consistently low by 5-8%. Therefore these were 
excluded, and the fits averaged over the remaining data. 
For February, it may therefore be that a similarly low 
(or high) result was obtained. Such a systematic error 
is nonetheless small compared to other uncertainties in 
our sensitivity calculation, so we do not carry this error 
through our calculation. 

2. Conversion to an effective V /in/ MHz threshold 

Using the piecewise linear approximations to k(v) it 
is possible to calculate the thresholds in V/m/MHz given 
the thresholds in a.d.u. These were constantly altered 
throughout the experiment. The thresholds varied signif- 
icantly between data channels and over time, since the 
thresholds were chosen to keep the trigger rate on each 
channel (rather than the thresholds themselves) constant. 

Since a positive detection requires a three-fold trigger, 



it is useful to define an effective signal detection thresh- 
old over all antennas. A good measure is to choose a 
mean V/m/MHz signal strength that over the bandwidth 
will have a 50% probability of triggering all three anten- 
nas. The random noise component can either increase 
or decrease the measured signal, so we require sufficient 
intrinsic signal strength above each individual antenna 
threshold that the chance of the random component push- 
ing the signal below threshold is small. Thus the effec- 
tive threshold is dependent upon this random component. 
The individual antenna thresholds were calculated from 
the mean recorded spectra from the relevant limb/center 
pointing calibrations by averaging the RMS signal over 
the entire 1 .024 GHz bandwidth, and converting the mea- 
sured RMS signal in a.d.u. into V/m/MHz using k. Dur- 
ing times of significant out-of-band RFI, the effective 
thresholds will vary due to a greater RMS signal, but 
since these occasions are both rare and have a low effec- 
tive efficiency, their contribution to the average threshold 
will be negligible and is neglected here. 

Assuming a normally-distributed RMS field strength, 
the probability of the total of signal plus noise falling 
above threshold for any given signal strength can 
be readily calculated, and thus the probability that 
the global condition (CA01A OR CA01B) and (CA03A OR 
CA03B) and (CA05A OR CA05B) will be met. The detec- 
tion probability is thus dependent on the alignment of the 
field vector with the A and B receivers - since the trig- 
ger condition is A OR B, the probability is highest when 
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FIG. 12: k(\) as measured for February 28''' for CA01 A based on lunar center/off-Moon, piece-wise fits to k(v). Dashed line 
shows the mean k(\). 



the field vector is parallel with either the A or B polar- 
ization directions, and lowest when it is 45° from both. 
Therefore the 'effective threshold' is defined for a sig- 
nal polarized at 22.5° (i.e. half-way between 0° and 45°) 
to either A or B. Since in general both the thresholds 
and RMS values are different for each polarization, we 
calculate the effective thresholds for signals polarized at 
both 22.5° (67.5°) and 67.5° (22.5°) to the A (B) receiver 
directions and average the results. By varying the raw 
signal strength until the calculated three-fold detection 
probability was 50%, the effective thresholds could be 
calculated. 



Doing so, we found that whereas each individual an- 
tenna trigger threshold was in the range 1.1-1.3 x 10~ 8 
V/m/MHz, the effective thresholds for a global trigger 
were in the range 1.45-1.6 x 10~ 8 V/m/MHz, i.e. an in- 
crease in threshold (decrease in sensitivity) of approx- 
imately 25%. In comparison, adding all three antenna 
voltages coherently (gain of V3 in signal-to-noise) and 
detecting at a higher level of Vthresh = 9.5Vrms (so the 
probability of a false detection for a ~ 30-hr experi- 
ment would be less than 0.1%; there would no longer 
be a coincidence check) would have produced an effec- 
tive threshold a few percentage less than the individual 
antenna thresholds, i.e. 1.05 — 1.25 V/m/MHz. Thus our 
inability to combine the signals coherently reduced our 
sensitivity by approximately 30%. 



3. Comparison with a possible experiment at Parkes 



An alternative instrument to the ATCA is the Parkes 
64 m single dish radio telescope, with an effective band- 
width of 300 MHz in the 1.2-1.5 GHz range. While 
the total sensitivity (here, area-bandwidth product) com- 
pared to the six ATCA antennas at 600 MHz bandwidth 
is 30% lower, with current technology we are unable to 
take advantage of the full ATCA collecting area, so we 
present a brief comparison of the two instruments. Ap- 
proximate scaling from ATCA to Parkes would suggest 
that: (i) we gain a factor of ^9 in area compared with 
one 22 m dish; (ii) we loose a factor of ^2 because we 
don't have a triple coincidence trigger and so must set the 
trigger threshold higher; (iii) we loose a factor of ~2 be- 
cause of the smaller bandwidth at Parkes; (iv) there is a 
modest gain because the lowest frequency is 1 GHz and 
not 1.2 GHz; (v) there is a modest gain because of the 
higher fraction of high quality RFI-free on-Moon time, 
though a lack of baseline may make RFI discrimination 
more difficult; (vi) there is a loss of a factor 2-3 be- 
cause the Parkes multibeam receiver can cover less of 
the lunar limb. Factors (i)-(iv) reduce the neutrino en- 
ergy threshold by ~ \/9/{2 x 2), i.e. overall the sensitiv- 
ity for our ongoing experiment at Parkes should be more 
than twice as sensitive as our 2008 experiment using the 
ATCA, while factors (v)-(vi) decrease the effective area 
to high-energy events by 50%. For targeted observations 
of potential sources of UHE particles (see Ref. []3~l|]), the 
Parkes telescope will be even more suitable, since for a 
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targeted observation, only part of the lunar limb would 
need to be observed. 



VII. RESULTS 
A. Search for lunar pulses 

The main search criteria used for eliminating false 
events were the timing requirements. The search win- 
dow is given by the apparent angular width of the Moon 
in the East- West direction, since the North-South com- 
ponent is unresolved by the East- West baseline. This is 
~ 0.5° at transit (when the Moon achieves its greatest el- 
evation), and considerably less near Moon-rise/set. This 
gives an intrinsic time-window of up to 23 ns (46 sam- 
ples) over the maximum baseline of 765 m. While nei- 
ther the raw nor the correlation-corrected times can be 
completely trusted, true lunar pulses will have a sharp 
time structure which will allow only a small variation 
in trigger times between antennas. Alternatively, those 
with extended structure (i.e. multi-peaked electromag- 
netic showers viewed away from the Cherenkov angle) 
will be due to the highest energy showers and therefore 
strong enough to give a correct correlation. 

Performing the search over both observation periods 
resulted in 60 candidates. Note that for any given lu- 
nar position there will always be two points on the hori- 
zon (one North and one South) which will have the same 
timing solution for an East- West baseline. Our search 
criterion is illustrated graphically in Fig. [T(j (see Sec. Q), 
where both the apparent direction of the Moon and the 
times and origins of candidate signals are plotted. Can- 
didate origins plotted for May 1 7 do not appear to be con- 
sistent with the lunar position because we have allowed a 
larger time window in the search due to an uncertainty in 
the timing calibration for that night (as previously noted). 
Note that all candidate events were detected during peri- 
ods of intense RFI. 

The candidate events were then searched through by 
eye for pulse-like events, and it was found that a majority 
of the candidates had a narrow-band RFI signature, with 
the recorded time-domain signals being strong over the 
entire buffer length. We did not use a more quantifiable 
measure than 'pulse-like' simply because the narrow- 
band RFI was so obvious and strong when present, and 
some 'pulses' had duration up to 30 ns. Since at these ex- 
tremely high energies, multiple cascade signatures (e.g. 
from nearby hadronic and purely electromagnetic cas- 
cades) might cause a lengthening of the expected pulse 
profile, we preferred to use timing criteria only if possi- 
ble. 

A minority of events - 16 in total - had a narrow time- 
structure, an example of which is given in Fig. All 
came within a two-hour period on February 27 th , which 
was one of the most RFI-intense periods of all the ob- 



servations. These could not be immediately excluded by 
eye, and had to pass more stringent tests. These are de- 
scribed below. 



1. Ensuring possible origin to within sampling accuracy. 

The search algorithm allowed both small deviations 
from a far-field event, and small offsets in direction from 
the Moon, to account for potential errors in the auto- 
mated alignment process. A check by eye of the cor- 
rected alignment times, if necessary including further ad- 
justments, would be expected to yield accurate timing 
information in cases where the detected event has signif- 
icant time-structure, as is the case with all 16 candidates. 
The correlation-corrected times occasionally needed fur- 
ther adjustment by one sample on a single baseline. For 
each candidate, the resulting alignment was compared vi- 
sually to that required for the event to have a far-field ori- 
gin; this was done quantitatively by comparing the buffer 
trigger times for CA03, £3, to that expected (£3) from £5 
and 1 1 : 

t' 3 = ((t 5 -h)(b 13 /b 15 )+h). (23) 

where b\ 3 and £15 are baselines between antennas CA01 
and CA03, and CA01 and CA05, respectively. An ex- 
ample of this procedure is given in Fig. [T4J. In no case 
did any of these events appear to be a far-field event. In 
most cases, structure was evident in both polarizations, 
but was so weak in one that the other only could be used 
for determining the alignment - however, an alignment 
to within sampling accuracy could always be obtained. 
Since at these times the Moon made an angle of nearly 
90° to the ATCA baseline, events within 360 km would 
result in a wavefront curvature with measurable differ- 
ences between £3 and £3 - in the case shown in Fig. |l4[ 
the distance is of order 100 km and so this event can not 
have a lunar origin. 

In principle we could have also used the expected dis- 
persion measure as a test to verify a pulse being of lunar 
or terrestrial origin. However, we deliberately chose ob- 
serving times when the dispersion was low, and so this 
reduced the power of such a test. Since our timing crite- 
ria were already sufficient to exclude all candidate events 
as RFI we did not pursue this approach. 



B. Effective apertures to an isotropic flux 

The simulation program described by James & 
Protheroe [|o[ was modified to weight the frequency 
spectrum output over the bandwidth by k(y)/k using the 
piece-wise linear approximation. The sensitivity of the 
experiment did not remain constant, but this was taken 
into account by running simulations using the lowest and 
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FIG. 13: A narrow time structure event from February 27 th . Note that the signal has been displaced vertically by ±30 as indicated 
for display purposes. 



highest values of thresh for each of the observation peri- 
ods. 

The deviation of the true TEC compared to that de- 
signed into our de-dispersion filters changed our sensi- 
tivity continuously. The average peak recorded signal 
strength as a fraction of intrinsic peak strength (Fig. [^) 
for each of the observation periods was folded into the 
simulations for different combinations of the two ex- 
treme Cherenkov spectra and the lowest, average, and 
highest sensitivity. 

It was found that changing the angles of the polarized 
receivers with respect to the lunar limb varied the calcu- 
lated isotropic apertures to neutrinos above 10 21 eV by 
less than 1%, and at most 4% at 10 20 eV, where the aper- 
tures are in any case very small. 

The two major uncertainties in this calculation are the 
UHE neutrino interaction cross-sections, and the effects 
of small-scale surface roughness (SSR), both of which 
are dealt with in detail in the appendices. In the follow- 
ing results, we do not incorporate the uncertainty in the 
cross-section for the simple reason that the uncertainty 
itself is so uncertain, since the scope for new physics 



at such high center-of-mass-energy collisions is large. 
Rather, we summarise the results of Appendix ^ by stat- 
ing that a doubling (halving) of the cross-section results 
in an increase (decrease) in the effective aperture by a 
factor of approximately 1.88. 

For the effects of small-scale surface roughness, we 
develop a toy model, and show both standard estimates 
("ATCA no SSR") and estimates adjusted for the results 
of calculations using this model ("ATCA SSR"). This 
model, while useful for understanding SSR phenomenol- 
ogy, is much less sophisticated than this topic requires, 
and the results based on it should be interpreted more 
as guides to indicate the nature of small-scale surface 
roughness effects. The important result is that whereas 
traditionally the effects of SSR were thought of as en- 
tirely negative, since the effect is to reduce the coherence 
of the wavefront over the cascade length, the resulting in- 
crease in angular spreading of the radiation makes a de- 
tection at the highest energies more likely. Therefore, the 
calculations including SSR have a lower effective aper- 
ture at low neutrino energies, and a much higher aperture 
at the highest energies. Future work is expected to pro- 
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FIG. 14: Comparison of best far field-fitted alignment (top) with the best unrestricted fit (bottom) for polarization B on one of the 
sixteen narrow-time-structure candidates on February 27 th , for the values of the offsets ( Afoi , A?5i, At^) shown. The 'wavefront 
curvature' is given by subtracting the measured ?3 from the expected value t £ given by Eq. E3l 



vide a more quantitative estimate of SSR effects. 



1. Apertures and limits to an isotropic flux 

The resulting range of effective apertures to an 
isotropic flux from each period is given in Fig. |l5[ 
assuming the presence of a radio-transparent mega- 
regolith. Also plotted are the effective apertures from 
prior experiments as calculated by James & Protheroe 
pc[]. For our ATCA observations, the threshold is 
lower than in past lunar Cherenkov experiments (since 
our high bandwidth has compensated for the smaller 
dishes), while the effective aperture at high energies is 
greater (due to increased coverage of the lunar limb and 
lower frequencies). Unlike the previous experiments 
with larger dishes, the sensitivity to UHE neutrinos is 
higher in the center-pointing mode (Feb 2008) than in the 
limb-pointing mode (May 2008), since the beam-width 
of ATCA near 1-2 GHz is comparable to the apparent 
diameter of the Moon. Note also that previous calcula- 
tions (Refs. [|5, 28 44]]) have not included the loss from 
a non-infinite sampling rate, and so their effective aper- 
ture should be reduced somewhat. 

The limit on an isotropic flux of neutrinos arising from 
our combined 2008 observations to a UHE v flux is given 
in Fig. |H as the band labeled "ATCA no SSR". Also 
shown is the limit - 'ATCA SSR" - when our toy model 
of small-scale surface roughness (see Appendix ph is in- 
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FIG. 15: (color online) The range of effective apertures (see 
text) for the LUNASKA February 2008 center-pointing (C) and 
May 2008 limb-pointing (L) ATCA observations, compared 
to that from previous experiments at Parkes ('Pks'), Kalyazin 
('Kal.'), and Goldstone ('GLUE'), in both limb- ('LA half- 
limb- ('HL'), and center- ('C') pointing modes (see PQ]), as- 
suming the existence of a sub-regolith layer of comparable di- 
electric properties to the regolith itself. We also show the effect 
on the aperture for the (center-pointing) Feb. 2008 observations 
including our adjusted toy model of small-scale surface rough- 
ness (C, SSR incl.) (see Appendix |b]), which is the reason for 
the abrupt (and artificial) change in aperture near 2 x 10 21 eV. 
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eluded, and the limits from previous experiments, includ- 
ing GLUE (dashed line labeled "GLUE") which is 
now believed to be approximately an order of magnitude 
too low as pointed out by James & Protheroe [ p0| ] and 
confirmed by Gayley et al. [^9| . The GLUE limit as re- 
vised upward by James & Protheroe is shown by the band 
line labeled "GLUE (JP09)". 



VIII. SUMMARY AND CONCLUSION 

In 2008 we carried out observations of the Moon at 
the Australia Telescope Compact Array using the lunar 
Cherenkov technique to search for the signatures of UHE 
neutrinos. Although no lunar pulses of any source were 
positively identified, the instantaneous apertures for all 
the observations were the most sensitive yet achieved 
using the technique. The corresponding limits on an 
isotropic flux are not as strong as those from RICE and 
ANITA, but our methods to improve sensitivity to certain 
patches of the sky - which, along with our limit on the 
UHE neutrino flux from Centaurus A and the Galactic 
Center [[34]], will be published separately - were a suc- 
cess: we found the exposure from our observations to 
Centaurus A to be higher than all previous experiments 
at neutrino energies of 10 22 -10 23 eV and above. 

The lack of detection of nanosecond lunar pulses re- 
ported here is consistent with limits set by the ANITA ex- 
periment. Importantly, our methods to discriminate be- 
tween true lunar pulses and RFI were very successful: we 
can be extremely confidant that no true lunar Cherenkov 
pulses were detected simultaneously by all three anten- 
nas, despite having of order ten million trigger events. 
This demonstrates primarily the power of nanosecond 
timing over a significant baseline, which was made pos- 
sible by our use of a very wide (here, 600 MHz) band- 
width, and was the main criterion used to discriminate 
against false events. Importantly, it was found in Sec. 
VII A 1 that for events of reasonable time -structure (the 



category into which all candidates must fall), the auto- 
matic procedures produced times accurate to better than 
1 ns. Therefore, relying on such a procedure to search for 
candidate pulses in the future is justified, making a com- 
parable analysis for a long observation run of a month or 
more feasible. 

Our first estimates of the effects of small-scale sur- 
face roughness on the detectability of lunar Cherenkov 
pulses at high frequencies — and the potential two orders 
of magnitude improvement in the effective aperture at 
the highest energies, and order-of-magnitude worsening 
of the threshold — only serves to emphasize the impor- 
tance of further work. Since these effects are indicative 
of 'messy' signals caused by roughness-induced interfer- 
ence along the length of the cascade, it also suggests that 
timing criteria on signal origin (i.e. that the signal comes 
from the Moon) be used instead of criteria on the pulse 



shape (i.e. that it is very close to impulsive). Our use of 
such criteria to eliminate all our candidate events proves 
that this is possible. It is important to note that the energy 
range at which small-scale surface-roughness effects are 
expected to increase the strength of our limit is the range 
at which our ATCA experiment is most competitive. 

For our observations, the use of an analog de- 
dispersion filter proved highly successful. The dispersion 
measure assumed in the filters' construction turned out 
to be very close to the actual values during observation 
periods, so that incorrect de-dispersion lost us less sensi- 
tivity than triggering inefficiency due to our non-infinite 
sampling rate. Though such filters must inevitably be 
superseded by a digital method, their continued use in 
the meantime will be valuable. Conversely, the finding 
that the loss from a finite sampling rate was greater for 
our simple trigger algorithm than that from incorrect de- 
dispersion is extremely important, and that in fact our 
'over' -sampling was an important factor in increasing 
(or rather, reducing the loss of) sensitivity. In future 
observations therefore the use of a smarter trigger algo- 
rithm that uses the already fully available information to 
reconstruct intermediate trigger values should be used, 
and perhaps should take as high a priority as digital de- 
dispersion. 

Compared with an alternative single-dish experiment 
at Parkes, our experiment at the ATCA provided more 
effective area to high neutrino energies at the expense 
of less sensitivity to lower energy events. Since the pa- 
rameter space at which the ATCA experiment is superior 
is best explored by low-frequency experiments such as 
NuMoon, we have transferred our efforts to utilising the 
Parkes dish, the results of which will be reported in a 
future contribution. We point out, however, that multi- 
telescope systems (such as the SKA and its pathfind- 
ers) will be more sensitive than single-dish experiments, 
and that our ATCA experiments would have been signif- 
icantly more sensitive had we not been limited to using 
only three antennas due to telescope upgrade delays. 

For future experiments at the ATCA or at other ra- 
dio telescope arrays, further improvements such as real- 
time coincidence logic between three or more antennas, 
or even the ultimate goal of a coherent addition of the 
signals, would also improve the sensitivity. Without a 
further analysis of the typical RFI structure, it is not pos- 
sible to determine the utility of real-time anti-RFI logic, 
though given the prevalence of RFI-triggered pulses, this 
too should be considered. 

The lessons learned above should in all cases be appli- 
cable to any use of the lunar Cherenkov technique with 
an array of radio antennas. The advantage of using a gi- 
ant radio array such as the SKA to search for lunar pulses 
has only been highlighted by these observations, espe- 
cially since it will be placed in a low-RFI environment. 

We have demonstrated techniques being developed by 
us ultimately for use with the Square Kilometer Ar- 
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FIG. 16: (color online) Model-independent 90% confidence limits, i.e. 2.3/[f e ffA e ff(£"v)] for our nominal 33.5 hr observations 
(effective on-time f e g- = 26.15 hr), on a total flux of UHE neutrinos (adjusted for all neutrino flavours) from our 2008 ATCA 
observations assuming a sub-regolith layer, both using the standard modeling ("ATCA no SSR") and using our adjusted toy model 
of small-scale surface roughness ("ATCA SSR incl."). In the 'SSR' case, the abrupt transition near 2 x 10 12 GeV is a model artefact 
(see Appendix §). Also, from previous experiments: GLUE [|l]]; IceCube Q; RICE Q]; ANITA FORTE NuMoon 
|p4]]; revised estimates by James & Protheroe [ pc| ] for Parkes, GLUE, and Kalyazin are shown by hatched bands (upper boundary 
- limit for 10 m regolith; lower boundary - 10 m regolith plus 2 km sub-regolith); Auger surface detectors [pi]]. The range on the 
band labeled ATCA reflects experimental uncertainties, while for previous experiments (where applicable) reflects the inclusion or 
otherwise of a sub-regolith layer. 



ray, and have been able with only 6 nights of obser- 
vations using the ATCA to produce the lowest limit to 
an isotropic UHE neutrino flux below 3 x 10 22 eV of 
any lunar Cherenkov experiment. While at present the 
isotropic limits from lunar Cherenkov experiments are 
not competitive with RICE Q, ANITA [f$| and (above 
3 x 10 22 eV) NuMoon, use of the SKA in several years' 
time for lunar Cherenkov observations will provide a 
very powerful technique for UHE neutrino astronomy 
[p0|] . With an estimated sensitivity to neutrinos 100 times 
less energetic than those detectable with our experiment 
at the ATCA, the SKA will be able to probe the as-yet 
untested predictions for a flux of neutrinos from the GZK 
process. Furthermore, our current experiment has been 
able to access regions of the sky not accessible to ANITA 
and NuMoon, and with better exposure than RICE above 
3 x 10 22 eV. Our flux limits to UHE neutrinos from Cen- 
taurus A will be discussed in a future paper. 



Appendix A: Variation of aperture with neutrino 
cross-section 



The greatest unknown in the calculation of the aper- 
ture is that of the UHE neutrino-nucleon cross-section 
Gviv, since this requires the extrapolation of experimental 
data over many orders of magnitude - hence the finding 
of Gandhi et al. [£2| that the cross-section could vary by 
a factor of 2 ± 1 at Kr° e V [ 52 1 . This conclusion was born 



out after a more recent calculation by Cooper-Sarkar and 
Sarkar [|53p using updated data from HERA estimated 
a neutrino-nucleon charged-current (CC) cross-section 
which is approximately 30% lower at 10 20 eV, and get- 
ting relatively smaller with increasing energy. There is 
also scope for 'new physics' to alter the cross-section 
even further. 

Given the range of uncertainty even within 'stan- 
dard physics', the determination of the UHE neutrino- 
nucleon cross-section is necessarily a scientific goal of 
UHE neutrino-detection experiments which is insepara- 
ble from the measurement of the UHE neutrino flux it- 
self, and in this context, instead of limits on the flux hav- 
ing a cross-sectional uncertainty, the limits should be set 
in flux-cross-section space. However, due both to con- 
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vention and the complexity of doing so, they are not. In- 
stead, to estimate the effect of uncertainty in the cross- 
section on our calculated effective area and flux limits, 
we calculate the fractional rate of change in the effec- 
tive aperture A e ff with the fractional rate of change in 



cross-section, i.e. 



In the limit where neutri- 



A,.ff cfoviv ' 

nos cannot penetrate a large part of the Moon and are 
seen only when they interact almost immediately in a thin 
layer at the Moon's surface ('down-going' interactions), 
we expect x^d^^J = 1' i- e - doubling the cross-section 
doubles the interaction rate, and vice-versa (a similar ef- 
fect is reached if the entire Moon is transparent to neu- 
trinos, but this is far from reality). Since this limit is 
approached for high neutrino energies observed at high 
frequencies, we expect results to be close to 1 . Any con- 
tribution from 'up-coming' neutrinos (neutrinos having 
penetrated a significant part of the Moon before reaching 
the surface) would reduce the result below 1, while no 
mechanism exists to increase the result above 1. 

To determine the sensitivity to the cross-section, we 
varied the cross-section by ±20%, and ran simulations 
for three primary neutrino energies - 10 21 , 10 22 , and 
10 23 eV - using the full range of ATCA sensitivities 
(best- and worst-cases) and configurations (February and 
May observations). We found the effect of varying the 
cross-section at a given energy was the same for all cases, 
and that over this range of a, was a l so constant 

at a given energy. Thus in Table |Vl|we give one value 
only of j 1 ^- for each primary neutrino energy, aver- 
aged over all cross-sections and observer configurations. 
That the values are very close to one confirms that down- 
going neutrinos dominate the detected interactions, espe- 
cially at high energies, so that a reduced (increased) es- 
timate for the neutrino-nucleon cross-section would pro- 
portionately decrease (increase) the ATCA experimen- 
tal aperture, with a corresponding increase (decrease) 
in the flux limit set from this experiment. Whether the 
cross-section can be deconvolved from the flux using 
(for instance) the average origin of the signal is a sub- 
ject for a future contribution. We also give in Table 
VI the reduced values of the charged-current neutrino- 
nucleon cross-section O v ncc using the fit given in Eq. 
3.5 of Cooper-Sarkar and Sarkar (CSS) and Gandhi et 
al. (GQRS), together with the implied reduction in ef- 
fective area A^ s /A^® RS under the assumption that the 
neutral-current cross-section scales with the charged- 
current cross-section in the CSS calculation. 



TABLE VI: The effects of changing cross-section on the sim- 
ulated effective aperture, assuming the neutral-current cross- 
section scales with the charged-current cross-section for the 
calculation of Cooper-Sarkar & Sarkar ySM. 



A c!! da vN 



J vNCC' "vWCC 

4 CSS IaCQRS 
^eff l A eS 



10 21 eV 


10 22 eV 


10 23 eV 


0.9 


0.915 


0.93 


0.55 


0.42 


0.30 


0.60 


0.46 


0.35 



since it takes the typical deviation at scales of order a 
wavelength, and treats it as a 'large-scale' (greater than 
a cascade length) phenomena. Here we define 'small- 
scale' roughness to be that on a scale between a cascade 
length (typically a few metres for hadronic cascades) and 
the wavelength in the medium. The consequences of re- 
stricting scales to this range is discussed later. The effect 
of small-scale roughness is to reduce the coherence be- 
tween radiation from different parts of the cascade. This 
in turn broadens the angular width over which the radia- 
tion is emitted while reducing the peak strength, and also 
allows transmission for angles of incidence greater than 
the critical angle, which is an effect observed for rough 
optical surfaces [p4[|. 

We carry out simulations using a toy model for an ex- 
treme case in which radiation from different parts of a 
shower in a near-surface cascade will see different sur- 
faces and thus be refracted semi-independently. We refer 
to the resulting effective aperture including our small- 
scale surface roughness model as A R S , while the effective 
aperture for the "standard case" without small-scale sur- 
face roughness is A^ ff . Because of interference effects, 
the true effective aperture is likely to be between these 
two extremes, and it may be reasonable to approximate 
this over a restricted frequency range by an "adjusted" 
effective aperture A^ ff = rA^ s + (1 — r)A R s . First we de- 
scribe the toy model, then the simulation method and the 
resulting aperture A R S . Interference effects and a method 
of determining r and hence A^ ff are dealt with, and fi- 
nally we discuss the validity of the approximations made 
and the necessity of a rigorous treatment of the effect of 
small-scale roughness. 



Appendix B: Variation of aperture with roughness model 

Current models of lunar surface roughness use a sin- 
gle surface slope over the entire length of the cascade 
through which radiation refracts. This is only an approx- 
imate treatment, as discussed by James & Protheroe J3(i|], 



1. Toy model calculations and results 

Assuming the smallest roughness scale to which ra- 
diation is sensitive is that of a wavelength, a cascade of 
length L might see up to Ng (~ L$/X) refractive surface 
elements, where L$ is the cascade length. Though it is 
unclear whether the applicable wavelength X is that of 
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the low-index or high-index medium, since our goal is 
to put an inclusive bound on the effects of such rough- 
ness, we choose X = X n , where A,,, is the wavelength in 
the medium of highest refractive index. We also choose X 
corresponding to the frequency / = (/ m in x /max) ' 5 - For 
our experiment at ATCA, / = 1 .47 GHz, so X„ = 11.8 cm 
in the regolith (n = 1 .73). 

We calculate Ns by scaling the energy- and medium- 
dependent L$ from its value of 12%o/p (12 radiation 
lengths given the medium density: ~ 4.7 m) for hadronic 
cascades in ice at 10 EeV fl55| ] by 7.5% per decade in en- 
ergy as per Williams [p4[], Since in our treatment, each 
portion of the cascade would 'see' a different refractive 
element, we break the cascade into Ns separate segments, 
rounding Ns up to take an integer value. Thus we arrive 
at a shower length given by Eq. Bl, and the number of 



shower segments Ns given by Eq. B2 



Ls 
N s 



12(Xb/p) 

(nf 
Ls 1 



1+0.075 log 



10 



E s 



10 EeV 



(Bl) 
(B2) 



We assume that each of the Ns segments contains an 
equal portion of the total excess tracklength of each cas- 
cade, so that the peak electric field amplitude from each 
is l/Ns that of the cascade as a whole. We also as- 
sume that the segment has length Ls/Ns, so the width 
of the Cherenkov cone is corr espondingly broadened by 
the same factor Ns. Table VII summarises the scaling re- 



lationships and number of shower segments for hadronic 
cascades at energies of 10 21 , 10 22 , and 10 23 eV in the 
regolith and megaregolith. 

The radiation from each cascade segment is treated in- 
dependently, with the simulated experiment able to de- 
tect none, all, or some of the cascade segments, with the 
detection probability of the primary neutrino being equal 
to one minus the probability that none of the individual 
cascade segments are detected. We calculate the trans- 
mission coefficients separately for each piece of surface, 
which results in the observed signal appearing to come 
mostly from those parts of the surface pointing roughly 
towards the observer. This method is identical to that 
described by James and Protheroe for handling cas- 
cades from interactions of secondary /j, T, and v. How- 
ever, it ignores the possibility of interference between 
the radiation from cascade segments, an effect which 
is unimportant for secondary interactions separated by 
large distances. We make an approximate adjustment for 
such interference b y ca lculating an 'adjusted' effective 
aperture A^ ff in Sec. B 2 - suffice to say for now that A^ ff 



must lie between A^ ff and A^ n . 

To generate a surface roughness deviate, we first gen- 
erate a deviate on the length scale Ls of the entire cascade 
as per James & Protheroe ftidj, where the slope tangents 
in each direction are sampled from a normal distribution 



with mean and standard deviation tanSnns given ac- 
cording to Eq. [B3] (modified from Olhoeft & Strangway 
by substituting Z^/lm for A,/lcm): 



tan S rl 



0.105L 



-0.22 



(B3) 



The distribution of surface roughness experienced by ra- 
diation at wavelength X„ (m) is found be replacing Ls by 
X n in the above equation. The extra roughness at small 
length scales can be thought of as additional small-scale 
deviation AtanS n ns superimposed over the large-scale 
roughness, which must produce the correct total value of 
S„n S (X n ). Since the small-scale deviation in this model is 
independent of the large-scale deviation, we find: 

tan 2 5 rms (A,„) = tan 2 ^^) + (Atan^s) 2 (B4) 



Substituting the formula for tanS'nns from Eq. B3 into Eq. 
B4 then gives the value of the additional surface deviate: 



Alans™ = 0.105 (X~ 



0.44 



-0.44 



0.5 



(B5) 



This gives a typical value for Atan.Srms of 8°; we approx- 
imate A tan S for each individual sub-showers to be in- 
dependent of each other for simplicity. Thus, for each 
cascade, we first calculate the length and generate a bulk 
surface normal by randomly sampling two slope tangents 
according to Eq. B3 , then for each sub-shower we modify 
the bulk surface normal by adding extra surface tangents 



sampled according to Eq. B5 



The additional complexity caused by breaking the cas- 
cade into segments means that we model only a simpli- 
fied version of our experiment, using three antennas with 
a fiat bandpass and single, circularly-polarised receivers 
operating in coincidence, with thresholds only approxi- 
mately that of our real experiment. We run the simulation 
both with (treatment described above) and without (stan- 
dard treatment) small-scale surface roughness for pri- 
mary neutrino energies in the range 10 20 <E V < 10 23 eV. 
Though the absolute values of the effective apertures A^ ff 
and Ag ff have little meaning for this fictional experiment, 
they are still illustrative to plot, which we do in Fig. |17| , 

Comparing A^ ff with A^ ff , the effect of small-scale 
roughness is significant. Since the peak emission from 
each cascade segment is a factor of Ns (~ 10) lower 
than that of the whole cascade, the effective neutrino 
energy detection threshold has been increased by a fac- 
tor of the same order. However, at the highest energies 
(> 10 22 eV), the probability of detection has increased 
by a factor of order 100, since the emission from each 
cascade is broader, and there are more cascades. In all, 
the effect of small-scale roughness on the detection prob- 
ability mimics that of observing at a lower frequency in 
the case when no such roughness is considered, though 
the expected time-domain signature would be quite dif- 
ferent. 
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TABLE VII: Parameters relating to the number of sub-showers in the surface roughness estimate, and some typical parameter 
values. Energies relate to the total hadronic cascade energy, which is typically 20% that of the primary neutrino energy. 



Layer 


Xo p X /p 
(g/cm 2 ) (g/cm 3 ) (cm) 


L S (m) 
10 21 eV 10 22 eV 10 23 eV 


Xj/n 
(m) 


N s 

10 21 eV 10 22 eV 10 23 eV 


Regolith 
Megaregolith 


22.59 1.8 12.55 
22.59 3.0 7.53 


1.51 1.62 1.73 
0.895 0.962 1.03 
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FIG. 17: Effective apertures of a fictitious experiment, calculated using large-scale roughness only (A^ ff ), our small-scale roughness 
approximation (Ag ff ), and an adjusted aperture which is a linear combination of them both (A^ ff - see Sec. B 2 ). 



2. Estimate of interference effects 

The aperture A* ff calculated by modeling small-scale 
surface roughness, as previously noted, excludes inter- 
ference between different cascade segments. Unlike ra- 
diation from two separate cascades, which if exiting the 
Moon in the same direction would be seen as two in- 
dependent signals by a detector, the arrival times of ra- 
diation from two cascade segments will likely be sep- 
arated by less than their duration - that is, they inter- 
fere. At one extreme ('case 1'), all the cascade segments 
will see the same surface and their radiation will exit in 
the same direction. Thus they interfere according to the 
standard Cherenkov condition, and the standard model- 
ing producing A^ ff is the correct treatment. At another 
extreme ('case 2'), the refracted radiation patterns from 
each cascade segment will not overlap, the cascades can 
be treated independently, and the results generated from 



our small-scale roughness model (i.e. A^ ff ) will be correct 
insofar as the surface generated is appropriate. 

We model our results as a linear combination of 
these two extremes by comparing the calculated aper- 
ture Agg^v) to that of A^ n (NsE v ). On average, the peak 
strength of the emission of each cascade in the A^. S (E V ) 
calculation will be the same as that of each cascade seg- 
ment in the A^ s (NsE v ) calculation, while the width of 
the radiation patterns from the segments under rough 
modeling will be k^Ns times greater, where kt allows 
for the slow growth of cascade length with energy: 

k L = (l+O^log!,,^)- 1 . (B6) 

Assuming no interference, the aperture in the rough 
'R' case will be one factor of Ns times larger than 
the standard 'S' case due to there being Ns as many 
sub-cascades, and a further factor of k^Ns d ue to 
the increased width as described above. Thus we 
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expect that the ratio A^ s (NsE v )/Al s (E v ) should be 
kLNsGH{NsE v )/Of]{Ev), where Gh is the cross-section 
for interactions producing hadronic cascades. This as- 
sumes that the interaction rate is proportional to the 
cross-section, which is accurate to order 10% (see Ap- 
pendix ^). We also assume that all hadronic interactions 
are dominated by neutrino-nucleon interactions (also an 
order 10% approximation Q), so that o H °c E° 363 [|§], 
and thus: 

A« s (N s E v )/Al s (E v ) = k L Nj 36 \ (B7) 

If the interference is complete (radiation always exits in 
an identical direction), and we assume a detection prob- 
ability of 1 within some part of the Cherenkov cone and 
otherwise so that there is no gain in detection probabil- 
ity from seeing the emission from two cascade segments, 
then only one of the Ns cascades will be detected in the 
A^ e (NsE v ) simulation, and we should find the ratio: 

A* s (NsEv)/A s eS (Ey) = k L N\ M \ (B8) 

We thus define the fractional overlap r such that r = in- 
dicates no interference (A^ ff applies), and r = 1 complete 
interference (A^ s applies), i.e.: 

A^(N s E v )/A s eS (E v ) = (1 - r)k L Nj M3 + rk L N l s M \B9) 

Using N s = 13 ± 1, we plot A* s (N s E v )/A s ef{ (E v ) in Fig. 
HH and fit for r to the high-energy regime where r is con- 
stant. Note that in this model the number N s of surface 
pieces refers to the number upon which the radiation is 
incident, not the number an observer sees. We expect an 
increasing ratio at low energies where a larger fraction 
of the emitted radiation is partially detectable (detection 
probability being neither nor 1), since then there is a 
gain in calculated aperture from two cascade segments 
radiating in the same direction. Though an increasing 
value of r is not observed at low energies - likely the ef- 
fect is obscured by the large uncertainties - we nonethe- 
less use only the range E > 10 21 eV for the fit, finding 
r = 0.70 ±0.06. 

An adjusted aperture A^ ff can then be calculated by 
A^ ff (£ v ) = rA s eff (E v ) + (1 - r)Af ff (£ v ). This has added to 
Fig. |l7] for our fictitious experiment. By assuming sim- 
ilar behaviour for our real experiment with the ATCA, 
we can calculate an adjusted aperture from our standard 
aperture only: 

Ai B (E v ) = 0.1A s eff (E v ) + 03k L Np 63 Al s (E v /N s ). 

(B10) 

The result has already been given for both our experi- 
mental aperture and limit in Figs. |l5| and [l6| respectively, 
by using the aforementioned values of r and Ns- 



3. Assessment of accuracy 

The model for small-scale surface roughness pre- 
sented here is intended as a toy model which deals only 
with the most significant possible effects of small-scale 
roughness at high frequencies. For instance, the sudden 
change in slope around 2 x 10 21 eV in Fig. |l7| is purely 
an artefact of our superposition of two extreme models 
- there is no reason to expect that a rigorous calcula- 
tion would produce any energy-dependentresults sharper 
than the slow (logarithmic) increase in cascade length 
with particle energy. 

In cases where the surface slopes are positively cor- 
related, and/or when a cascade is sufficiently deep that 
the far-field conditions are satisfied at the surface over 
most of the cascade, then our model of treating the ra- 
diation from each cascade/surface piece independently 
breaks down. While our introduction of the overlap pa- 
rameter 'r' goes some way to correcting for such cor- 
relation, in these cases the value of 'r' would likely be 
larger, so the effects of small-scale roughness estimated 
by our toy model are more likely to be over-estimates 
than under-estimates. 

Our model also ignores roughness on scales smaller 
than a wavelength, which can not be dealt with by split- 
ting the cascade into segments as in our toy model (it 
is not possible to produce multiple-refraction effects for 
structures of size less than that of a wavelength). How- 
ever, we do expect such roughness to broaden the emis- 
sion at the expense of peak refracted field strength - for 
our current model, this would mostly cause an increase 
in the overlap ratio r in our results. Therefore, we ignore 
sub-wavelength roughness until a more complete model 
can be constructed, which we leave to a future work. An 
additional approximation is that we have assumed only 
a single frequency for our surface-roughness estimates, 
rather than a continuous range. 

Despite these shortcoming, we have developed the first 
treatment of roughness on scales smaller than a cascade 
length, an important effect for high-frequency observa- 
tions that has been ignored in all previous calculations. 
Whether this effect helps or hinders neutrino detection 
will depend upon the shape of the UHE neutrino spec- 
trum. What we can say is that the difference in detec- 
tion probability between high- and low-frequency obser- 
vations may not be as dramatic as previously thought, 
but that a reconstruction of cascade parameters from a 
detected event may prove more difficult. 
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FIG. 18: The ratio A^ tt (N s E v )/A s ef[ (E w ). The error bars reflect both uncertainties in the apertures themselves from the Monte Carlo 
simulation (inner error bars), and also including the error from Ns = 13 ± 1. The fit is to the range E v > 10 21 eV only. 
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